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Abstract 

"We investigate the imaginary parts of the nucleon electromagnetic and axial 
form factors close to threshold in the framework of heavy baryon chiral per- 
turbation theory. For the isovector electromagnetic form factors, we recover 
the well known strong enhancement near threshold. For the isoscalar ones, we 
show that there is no visible enhancement due to the three-pion continuum. 
This justifies the use of vector meson poles only in dispersion-theoretical cal- 
culations. We also calculate the imaginary part of the nucleon isovector axial 
form factor and show that it is small in the threshold region. 
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I. INTRODUCTION AND SUMMARY 



The electromagnetic and axial structure of the nucleon as revealed e.g. in elastic electron- 
nucleon and (anti)neutrino-nucleon scattering is parameterized in terms of the six form 
factors Fi^2{t) C^A,p(i)0 (with t the squared momentum transfer). The understanding of 
these form factors is of utmost importance in any theory or model of the strong interactions. 
Abundant data on these form factors over a large range of momentum transfer exist and this 
data base will increase when the recent/upcoming experiments performed at ELSA, MAMI 
and TJNAL will have been analyzed. Dispersion theory is a tool to interpret (and cross 
check) these data in a largely model-independent fashion. The form factors can be written in 
terms of unsubtracted dispersion relations and their absorptive parts are often parameterized 
in terms of a few vector meson poles. This procedure is based on the successful vector 
meson dominance (VMD) hypothesis, which states that a photon (or W/Z-hoson) couples 
to hadrons only via intermediate vector mesons. However, as already pointed out in 1959 
by Frazer and Fulco such an approach is not in conformity with general constraints from 
unitarity and analyticity. In particular the singularity structure of the triangle diagram is 
not respected at all when using only vector meson poles. As a consequence of this singularity 
structure the two-pion continuum has a pronounced effect on the isovector spectral functions 
on the left wing of the p-resonance. This becomes particularly visible in the determinations 
of the corresponding nucleon mean square radii. This effect was quantified by the Karlsruhe- 
Helsinki group in their seminal work on the nucleon electromagnetic form factors . This 
analysis was recently refined by the Mainz-Bonn group accounting for new data and the 
high energy constraints from perturbative QCD However, there still remains one open 
end concerning the dispersion-theoretical analysis. In the isoscalar electromagnetic channel, 
it is believed (but not proven) that the pertinent spectral functions rise smoothly from the 
three-pion threshold to the tu-meson peak, i.e. that there is no pronounced effect from the 
three-pion cut on the left wing of the cj-resonance (which also has a much smaller width 
than the p-meson). Chiral perturbation theory can be used to settle this issue. An 
investigation of the isoscalar spectral functions based on pion scattering data and dispersion 
theory as done for the isovector spectral functions PQ] seems not to be feasible at the 
moment since it requires the full dispersion-theoretical analysis of the three-body processes 
ttN — >■ ttttN (or of the data on NN Svr). 

The one loop calculation of the isovector nucleon form factors indeed shows the strong 
unitarity correction on the left wing of the p-meson (i.e. slightly above threshold) For 
the isoscalar form factors, a calculation of the imaginary parts of certain two-loop diagrams 
will reveal whether there is some enhancement above t = 9M^ or justify the common 
assumption that one has a smooth isoscalar spectral functions driven by the tu-meson at low 
t. Similarly, the spectral function of the isovector axial form factor is supposedly dominated 
by correlated vrp-exchange (or the ai-meson) in the region above threshold. Nothing is yet 
known about a possible strong enhancement of the pertinent imaginary part slightly above 
threshold. It is exactly this problem we wish to address in this paper by evaluating the 
imaginary parts of the pertinent two loop diagrams at chiral order and contributing to 
the isoscalar electromagnetic and the isovector axial form factors in the framework of heavy 
baryon chiral perturbation theory. 



'^■^In what follows, we will not be concerned with the induced pseudoscalar form factor Gp{t) 
which is dominated by its pion-pole term Ag^f^mF^/^M^ — t). 
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The pertinent results of this investigation can be summarized as follows: 



(i) We have calculated the imaginary parts of the isovector electromagnetic form factors in 
heavy baryon CHPT to order (in the effective Lagrangian), thus extending previous 
studies The strong enhancement slightly above threshold on the left wing of 

the p-resonance due to the two-pion continuum is recovered and good agreement with 
the empirical analysis is found if one subtracts from it the p-meson contribution (by 



setting the pion charge form factor F^(t) = 1) as described in section |T|). 



(ii) We have determined the leading contribution to the imaginary parts of the isoscalar 
electromagnetic form factors in the threshold region by evaluating the imaginary part of 
the relevant two-loop diagrams at order q' . As previously anticipated, these imaginary 
parts are very small and rise smoothly with increasing t. The common procedure 
[^] of describing the isoscalar electromagnetic spectral functions by the cj-meson pole 
(and higher mass poles and/or continua) is therefore justified. 

(iii) We have also calculated the imaginary part of the nucleon isovector axial form factor 
in the threshold region. It shows a behavior similar to the isoscalar electromagnetic 
spectral functions, i.e. a smooth rise from threshold. Furthermore, this three-pion 
contribution is numerically small. 



II. BASIC CONCEPTS 

In this section, we assemble all necessary definitions for the calculation of the imaginary 
parts (spectral distributions) discussed below. This section serves mainly the purpose of 
fixing the notation underlying our analysis. This material is mostly not new but necessary 
to keep the manuscript self-contained. The reader familiar with it is invited to skip this 
section. 

Consider first the nucleon electromagnetic form factors. They are defined by the nucleon 
matrix element of the quark electromagnetic current, 

{N{p')\q^^Qq\N{p))=u{p') 

with t = {p' — pY the invariant momentum transfer squared and Q the quark charge ma- 
trix, Q = diag(|, — |, — |). Fi{t) and F2{t) are called the Dirac and Pauli form factors, in 
order. Following the conventions of [Q, the electromagnetic form factors are decomposed 
into isoscalar (5) and isovector (V) components 

F,,,{t) = F,%{t) + T'Fr^,{t) , (2) 

subject to the normalization 

(0) = F^O) = \ , Fi{Q) = , F^iO) = , (3) 

with Kp = 1.793 [nn = —1.913) the anomalous magnetic moment of the proton (neutron). 
In what follows, we will work with the electric and magnetic (Sachs) form factors. 



rF,{t) + —a>^-'ip'~p),F2{t) 



u{p) 



(1) 
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G'/it) = Ff'^(t) + ^i^2^'^W , G'/it) = Ff'^(t) + (4) 

in the isospin basis, and m = 938.27 MeV denotes the nucleon mass. 

Similarly, the nucleon matrix element of the quark isovector axial current is parameter- 
ized in terms of two form factors. 



{N{p')\qYl,T'^q\N{p)) = u{p') YGA{t) + ^(p' - p)^ Gp(t) 



^^r'^uip) , (5) 



with Gjiit) called the axial and Gp{t) the induced pseudoscalar form factor. In what follows, 
we will only be concerned with the axial form factor. Its normalization is given by G'yi(O) = 
QA = 1.26. 

We briefly discuss now the pertinent dispersion relations. Let F{t) be a generic symbol 
for any one of the five nucleon electroweak form factors. We assume the validity of an 
unsubtracted dispersion relation of the form 

Fit) = - / ^ dt' , (6 

where the —ie in the numerator is necessary ii t > Iq, since in that region F{t) is complex 
valued. The spectral function Im F{t) is different from zero along the cut from to to +oo 
on the real t-axis, with to = 9M^ for the isoscalar electromagnetic and the axial case. 
M,r = 139.57 MeV denotes the pion mass and we neglect isospin breaking effects in the pion 
and nucleon masses. For the isovector electromagnetic form factors the threshold lies at 
to = 4M^. The proof of the validity of such dispersion relations in QCD has not yet been 
given [§]. Eq.(^ means that the electroweak structure of the nucleon is entirely determined 
from its absorptive behavior, with the data for F{t) given (mainly) for t < 0. In what 
follows, we will concentrate on the threshold behavior of the various spectral functions, i.e. 
the range < t < 30 M^. The upper end of this range corresponds to the mass of the 
first isoscalar and isovector vector mesons, i.e. the a;- and the p-meson (M^ = 31.4 M^, 
Ml = 30.3 M^). 
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III. ISOVECTOR ELECTROMAGNETIC SPECTRAL FUNCTIONS 

Unitarity determines the spectral functions of the isovector electromagnetic form factors 
in the interval 4M^ < t < 16M^ uniquely as,[ 



lu.Gl{t)=^^^^fl{t)Fr{t)^ 

lrnGUt)=^^^^^fAt)Fr{t), (7) 



with F^{t) the pion charge form factor and the /j-(t) are the P-wave ttN partial wave 
amplitudes in the t-channel (vrvr — >■ NN), cf. fig. 1. The partial wave amplitudes f±{t) have 



#^In practise, this representation holds up to t ~ 50 
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a logarithmic singularity on the second Riemann sheet (originating from the projection 
of the nucleon pole terms in the invariant nN scattering amplitudes) located at 

= - M^/m^ = 3.98 Ml , (8) 

very close to the physical threshold at to = 4M^. The isovector form factors inherit this 
singularity (on the second sheet) and the closeness to the physical threshold leads to the 
pronounced enhancement of the isovector spectral functions weighted with as shown 
in fig. 2. Note that the ratio of the magnetic to the electric spectral function is very close 
to 1 + 4.71, i.e. the ratio of the form factors in dipole approximation. For 

later comparison, we also show the corresponding curves with F^{t) = 1 (to suppress the 
contribution from the p-meson). This strong enhancement is obviously very important for 
a precise determination of the nucleon isovector mean square radii . 

In the framework of relativistic baryon CHPT, the isovector spectral functions were 
calculated in The correct analytic structure naturally emerged in this calculation and 
thus the strong threshold enhancement. However, in that framework there is no consistent 
power counting scheme due to the extra energy scale related to the nucleon mass in the chiral 
limit. This can be overcome in the heavy baryon approach, which is discussed in detail in 
the review 0. In that paper, the imaginary parts of the isovector Dirac and Pauli form 
factors were calculated to order and compared to the result of 0. Here, we go one step 
further and calculate also the first corrections to these results at order q^, thus completing 
the one loop calculation. The imaginary part of the isovector electric form factor can be 
given entirely in terms of lowest order parameters 



ImG^(t) 



t - AMI r 



t - AMI + Ai^i - 



2\ 



9\{t-2Ml 



2\2 



(9) 



where the term in the square brackets is the order result and = 92.4 MeV is the pion 
decay constant. We remark that in the expressions for the imaginary parts we make use of 
the Goldberger-Treiman relation to determine the value of qa-, Qa = g-KNFT^/fn = 1.32, since 
here the pion-nucleon coupling constant g.„N = 13.4 appears naturally and not the axial 
current coupling qa- The expression for the imaginary part of the isovector magnetic form 
factor to order contains the low-energy constant C4 of the dimension two pion-nucleon 
Lagrangian (we take the value C4 = 2.25 GeV~^ determined in [0] from some vrA^ data). 



lmGl,{t) 



V u^ 9\m{t - AMI) , " 



64F2v^ 1927rF2v^ 



;i + 4mc4)(t - AMI) + 9\{lQMl - 7t) 



(10) 



We remark that at this order the normal threshold to = AM^. and the anomalous threshold 
tc = 4M^ — M^/rn?' on the second Riemann sheet still coalesce (it requires accuracy q^ to 
separate to and tc). As a consequence of this, the isovector imaginary parts show an abnormal 
threshold behavior in the 1/m expansion, i.e. they do not start out as (t — 4M2)^/2 g^g jj-^ 
eq.(7). The coalescence of both thresholds leads to smaller exponents. This problem is, 
however, largely academic as the comparison of the so calculated imaginary parts (weighted 
with l/t^) shown in fig. 3 with the dashed-dotted lines of fig. 2 reveals (one should only 
compare to the empirical spectral functions with the contribution of the p-meson subtracted 
since this resonance is not included in the effective field theory). In fig. 3, the results of 
the order q^ and order calculations are shown separately by the dashed and solid lines, 
respectively, indicating that the empirical observation of the strong peak close to threshold 
can indeed be explained within the heavy mass expansion of the chiral effective pion-nucleon 
field theory. This shows that the method can be used to investigate the isoscalar spectral 
distributions in the threshold region (below the a;-resonance) . 
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IV. ISOSCALAR ELECTROMAGNETIC SPECTRAL FUNCTIONS 



The imaginary parts of the isoscalar electromagnetic form factors open at the three-pion 
threshold to = 9 M^. The three-pion cut contribution is depicted in fig. 4. In the NN center- 
of-mass system, the pertinent anti-nucleon and nucleon four-momenta are Pi = (-\/t/2,p) 
and P2 = {\/i/2, —p). The analytic continuation of the three-momentum p into the region 
t < 4m^ takes the form 



p 



p^ 



w? = imp + 



(11) 



in the heavy mass expansion, with p a real unit vector. Application of the Cutkosky rules to 
three-pion intermediate state with four-momenta li, I2 and h= P1+P2 — I1 — h (cf. fig. 4) 
leads to 



Im 



A 



i(27r)^ 



d% 



1 



1 



1 



z(2vr)4 Ml - II Ml - q Ml - q 



B 



ydVsiAB), (12) 



where the symbol 'A' refers to the 7 — Svr and 'B' to the Sir NN transition, respectively, 
and dr^ is the measure on the invariant three-body phase space. In what follows, we have 
to evaluate integrals over the three-pion phase space. Tensorial integrals can be reduced to 
scalar ones by the following formula (or any permutation of it in the indices '1' and '2') 



dV^H{...)rJ{il = dVsH 



X 



1 - x')5'^p'' + -{z - xy){5'^p> + 5^''p 



"ik ■ 



xz) p' p^ 



(13) 



with 



0J7 



Ml , , 
UJ1UJ2 



= 1,2) , x = k-p , y 

Vi{uj^ + UJ2) + Ut + Ml) 



k -P 



\k\ \k 



(14) 



and the function H{. . .) depends on the variables ui, U2, x, y, z, t. The integration over the 
three-pion phase space can be expressed as a four- dimensional integral of the form 



dT.H(...] 
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^ f f duiduJ2 I I 



dxdy 



E y/1 



y2 _ ~;2 _|_ 2xyz 



Hi...), (15) 



+ 2xyz > 0. 



where the interior of the ellipse E is defined by the equation 1 — x"^ — y^ 
The elhpse has semi-axes of length \/l — z and a/I + z, respectively. 

In chiral perturbation theory, to leading order g^, the two-loop diagrams shown in fig. 5 
can contribute to the isoscalar imaginary parts, but graph (d) vanishes because of an isospin 
factor zero. All (relativistic) Feynman rules needed to calculate these graphs can be recov- 
ered from appendix A of with the exception of the anomalous 787? vertex. In the chiral 
effective Lagrangian approach, it follows from the Wess-Zumino-Witten term (for a review, 
see m). 
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C^l = A, Tr (a.f/f/t d^UU^ dpUU^ ) , (16) 

with the photon field and U = 1 + ir • vf/F + . . . collects the pions. Note that this vertex 
has dimension four and therefore the two-loop graphs for the imaginary parts start only at 
order . The corresponding Feynman insertion (with the sign convention e"^^^ = —1) is 

737r- vertex: - e ■ [h x h) e'^^'^ (17) 

where 'a, 6, d are pion isospin indices and e is the photon polarization vector. With these 
rules at hand and using eqs. (11-15) one arrives after straightforward but tedious calculations 
at the isoscalar imaginary parts in the limit m ^ oo. We give only the results in this limit 
since these represent the genuine leading order contributions with all higher order effects 
(starting at order q^) switched off. Furthermore, only in this limit the angular integrations 
in eq.(|l5|) can be performed analytically. We find: 



ImCKt) = -|44^ / / duj,du2 \h\ U Vl^ arccos(-z) , (18) 



ImGiit) 
+W{t) 



gAin 



3f 



(8vr)4F6 

+g\{bt^ + m^/'^M^ - Ult^Ml + m^''^Ml + 277tM^ - bSVtM^ - 135M^) 



(19) 



with 



m 

W{t) 



1^ Vt-M^ + - 2ViM^ - 3M2 



2t3/2 

Vi-M, 
96^3/2 



t-2v^M^-3M2 . 



(20) 
(21) 



Note that in the infinite nucleon mass limit ImG^(t) comes solely from graph (c) in fig. 5 
and quite astonishingly one can evaluate all integrals in closed form for ImGf^(t). For the 
sake of completeness we give also the isoscalar electromagnetic spectral functions in the 
chiral limit = 0. The resulting expressions can be given in closed form, 



o s 



Im Ge (t) 



1 f 9At 
105 U67rF2 



{t) - 6(167r)4F6 



(22) 



and exhibit a simple power like dependence on the variable t. In eq. (p2D all quantities are 
to be taken at their values in the chiral limit, as denoted by the with the exception of 
the pion decay constant whose chiral limit value is called F. The behavior near threshold 
to = 9M2 of the imaginary parts for finite pion mass, eqs.(p!^,p!9|), is 
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ImG'l(t) ~ (yt-3M^)^ ImG'f^(t) ~ (Vt-3M^)^/^ (23) 
which corresponds to a stronger growth than pure phase space 

duJidiU2 \h X Tap ~ (v^ - 3M^)^ . (24) 



22<1 



This feature indicates (as in the isovector case) that in the heavy nucleon mass hmit m ^ oo 
normal and anomalous thresholds coincide. In order to find these singularities for finite 
nucleon mass m an investigation of Landau equations is necessary [|lT|. By using standard 



techniques |TTI] we are able find (at least) one anomalous threshold of diagrams (a) and (b) 

at 

= (^^A - M^/m^ + ^1 - M2/m2) , t, = 8.90M^ (25) 

which is very near to the (normal) threshold to = 9 and indeed coalesces with to in 
the infinite nucleon mass limit. We note that diagram (d) does not possess this anomalous 
threshold tc = 8.90 M^, but only the normal one. We do not want to go here deeper into the 
rather complicated analysis of the full singularity structure of all two-loop diagrams in fig. 5 
but are mainly interested in the magnitude of the isoscalar electromagnetic imaginary parts. 
The resulting spectral distributions again weighted with l/t^ are shown in fig. 6. They show 
a smooth rise and are two orders of magnitude smaller than the corresponding isovector ones, 
cf. fig. 3. This smallness justifies the procedure in the dispersion-theoretical analysis like 
in m to describe the isoscalar spectral functions solely by vector meson poles starting with 
the a;-meson in the low energy region. Nevertheless, it may be worthwhile to include these 
calculated isoscalar imaginary parts in future dispersion analyses. We finally remark, that 
Im G%i^j{t)/t^ which have the same asymptotic behavior (for t oo) as Im G^j^{t)/t^ 
(considering only the leading contribution) do still not show any strong peak below 
the cj-resonance. Im G^{t)/t'^ is monotonically increasing from to = 9M^ to t = 30 
and Im Gf^(t)/t'^ develops some plateau between t = 20 and 30 M^. This observation is 
a further indication that there is indeed no enhancement of the isoscalar electromagnetic 
spectral function near threshold. Even though the isoscalar and isovector electromagnetic 
form factors behave formally very similar concerning the existence of anomalous thresholds tc 
very close to the normal thresholds to, the infiuence of these on the physical spectral functions 
is rather different for the two cases. Only in the isovector case a strong enhancement is 
visible. This is presumably due to the different phase space factors, which are (t — to)^/^ 
and (t — to)^ for the isovector and isoscalar case, respectively. In latter case, the anomalous 
threshold at tc = 8.9 is thus effectively masked. 



V. ISOVECTOR AXIAL SPECTRAL FUNCTION 

The calculation of the imaginary part of the axial form factor Im G'a(^) proceeds along 
the same lines as outlined in the previous section. One has to consider the same graphs 
shown in fig. 5 (with the wiggly line denoting now the external axial source, i.e. the W- 
boson). All four graphs contribute and performing the necessary integrations, the imaginary 
part of G'a(^) in the heavy mass limit reads 
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— * 



arccos — 2; 



(26) 



In the chiral limit M^^ = the isovector axial spectral function shows a simple dependence 
of the form 



Im Ga (t) 



9a t' 



o2 



22 



(167r)3F4 

We remark again that the behavior near threshold, 

ImG'^(t) ~ {yft-Ud^f 
is a stronger rise than pure phase space 



K = 6.663 . 



'<i 



duJidu2 1/1 1 1/2 



t - 3M^ 



(27) 



(28) 



(29) 



If one considers individual diagrams one finds that only the graphs (a), (b) and (c) which have 
intermediate nucleon propagators show the stronger threshold behavior of eq.(^). Diagram 
(d) without a nucleon propagator follows the phase space argument eg. (|29|) . The obvious 
reason for these features is the coincidence of normal and anomalous threshold in graphs (a), 
(b) and (c) which does not occur in diagram (d) in the limit m — >■ 00 Finally, fig. 7 shows the 
isovector axial spectral function Im GA{t) which again has a smooth rise from threshold and 
is about one order of magnitude smaller than the isovector electric spectral function. Even 
though Im GA{t) starts at order it is only three times as large as the isoscalar magnetic 
spectral function lmGfj{t) of order g^. We remark that Im GA{t)/t^, which has the same 
asymptotics as Im G^{t)/t'^ (to order g^), is still monotonically increasing from threshold to 
t = 30M2, indicating once more that the three-pion continuum has no pronounced effect on 
the axial spectral function near threshold. Of course, one also has to keep in mind that the 
available data for Ga(^) are not yet that precise to perform accurate dispersion-theoretical 
calculation. 

This completes our analysis of the spectral functions related to the isoscalar electromag- 
netic and isovector axial nucleon form factors near threshold. We do not find a visible strong 
enhancement slightly above threshold = 9 as in the isovector electromagnetic case, 
even though these imaginary parts do not follow the phase space in the heavy nucleon mass 
limit. The presence of a nearby anomalous threshold at tc = 8.9 has no substantial effect 
on the strength of the three-pion continuum. 
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Figure captions 

Fig.l Two-pion cut contribution to the nucleon isovector electromagnetic form factors. Wig- 
gly, dashed and solid lines represent photons, pions and nucleons, in order. 

Fig. 2 Spectral distribution of the isovector electric and magnetic nucleon form factors 
weighted with l/t"^ on the left wing of the p-resonance as calculated in Shown 
are lmG\^(t)/t'^ (upper solid line) and ImG^(t)/t^ (lower solid line). Suppressing the 
p-meson contribution as described in the text {F^(t) = 1) leads to the respective 
dashed-dotted lines. 

Fig. 3 Spectral distribution of the isovector electric and magnetic nucleon form factors 
weighted with l/t^ calculated in heavy baryon CHPT. Shown are lmG\^(t)/t'^ (upper 
lines) and ImG^(t)/t^ (lower lines). The solid and dashed lines refer to the order 
and calculations, respectively. 

Fig. 4 Three-pion cut contribution to the nucleon isoscalar electromagnetic form factors. 
The vertex denoted A' describes the 7 ^ 37r coupling and the one denoted 'B' the 
377 NN transition. 

Fig. 5 Two-loop diagrams contributing to the imaginary parts of the isoscalar electromag- 
netic and isovector axial nucleon form factors. 

Fig. 6 Spectral distribution of the isoscalar electric and magnetic nucleon form factors 
weighted with 1/t^ in the heavy nucleon limit. Shown are ImG'f^(t)/t^ (upper line) 
and ImG|(t)/t^ (lower line). 

Fig. 7 Spectral distribution of the isovector axial nucleon form factor weighted with l/t"^ in 
heavy nucleon limit. 
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